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I. INTRODUCTION 

The phenomenon of entanglement jj], 0, 13 is a remarkable feature of quantum physics that has been identified as 
a key ingredient in many areas of quantum information theory including quantum key distribution Q, superdense 
coding [f| and teleportation [|| . However the general problem of how to quantify Q the level of entanglement in an 
arbitrary multipartite system remains unresolved. There has been some progress towards a solution |7l IsL Ifl HoL ITU 
121 Hil Ha Ha . Il7| . but the task at hand is generally considered a difficult one and may never be completed. We 
are thus lead to consider simple computable measures of entanglement 0, 0] that although cannot fully characterize 
the multipartite nature of the correlations, may nevertheless still provide a useful gauge of their levels. 

In this article we investigate the average bipartite entanglement, over all possible divisions of a multipartite system, 
as a useful measure of multipartite entanglement. Such measures might be considered the least sophisticated of choices; 
however, their simplicity allows theoretical calculations to be exercised with ease. We will restrict our study to pure- 
state entanglement where the subsystem linear entropy is a clear choice for the bipartite measure. It was recently 
shown by Brennen [23] that an entanglement measure proposed by Meyer and Wallach is of the above described 
form, and hence, the multipartite entanglement measures considered in this paper may be viewed as generalizations of 
the Meyer- Wallach measure. Our measures may also be viewed as variations of those considered by Pope and Milburn 
|2l| where instead the minimum bipartite entanglement was considered. 

We show that the average bipartite entanglement elects self-dual quantum-error-correcting codes to the status of 
maximally entangled states. The connection between entanglement and quantum-error-correcting codes has been 
highlighted elsewhere (e.g. H3); however we make this relationship explicit by expressing the average entanglement 
of encoded states in terms of the weight distribution of the code. We also investi gate the multipartite entangling 
power of quantum evolutions. A simple extension of the work of Zanardi et al. 23] allows the derivation of an 
explicit formula. Such formulae are relevant to current studies in the entangling capabilities of chaotic systems 
ISIillillillMliSlilEi An example treated in this article is the quantum kicked 

rotor. 

The paper is organized as follows. In the next section we introduce the Meyer- Wallach entanglement measure and 
its generalizations. The connection between these measures and quantum-error-correcting codes is discussed in Section 
IIIII This relationship is further strengthened in Section llVI where we derive a formula for the average entanglement 
over a subspace. In Section we derive a formula for the multipartite entangling power of an arbitrary unitary. 
Finally in Section IVII we conclude by applying our results to the quantum kicked rotor. 



II. A CLASS OF MULTIPARTITE ENTANGLEMENT MEASURES 



It is generally accepted that when a bipartite quantum system is in an overall pure state, there is an essentially 
unique resource-based measure of entanglement between the two subsystems. This measure is given by the von 
Neumann entropy of the marginal density operators To ease theoretical calculations, one often replaces the 
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von Neumann entropy with its linearized version, the linear entropy. For a bipartite system in an overall pure state 
\ip) £ C Da (8 <C Db , the subsystem linear entropy is defined as 

S L (iP)=rt(l-tTp A 2 ) p A =tt B \il>M (1) 

where the normalization factor, n = D/(D — 1) with D = min(D J 4, Db), is chosen such that < Sl < 1. The state 
is separable if and only if Si = 0, and maximally entangled when Sl — 1 • 

In general, as the number of subsystems increases, an exponential number of independent measures is needed to 
quantify fully the amount entanglement in a multipartite system. Consequently, the following entanglement measures 
cannot be thought of as unique. Different measures will capture different aspects of multipartite entanglement. 

The Meyer- Wallach measure |l9|, Q(ip), which can only be applied to multi-qubit pure states \ip) £ (C 2 )®", is 
defined as follows. For each j = 1, . . . ,n and b £ {0, 1}, we define the linear map ij(b) : (C 2 )®" — > (C 2 )® n_1 through 
its action on the product basis, 

ij(b)\xi) (8> • • • (8> \x n ) = S bxi \xi) • • • <g> \xj_t) <g> \x j+ i) ® • • • (g) \x n ) (2) 

where Xi £ {0, 1}. The Meyer- Wallach entanglement measure is then 

Q(i>) = - £r>(v(o)|v), (3) 
"i=i 

where 



D(^),|^)) = -|W>| 2 - (4) 

Meyer and Wallach showed that Q is invariant under local unitary transformations and that < Q < 1, with Q (^) = 
if and only if is a product state. 

Recently, it was shown by Brennen |2f)j | that Q is simply the average subsystem linear entropy of the constituent 
qubits: 

QW = 2^1-i^trp fe 2 ^ (5) 



where pk is the density operator for the fc-th qubit after tracing out the rest. This simplification is easily understood 
|42| by first showing that D(ij(0)\ip), ij is unchanged by a local unitary applied to the j-th qubit (a fact already 
proven by Meyer and Wallach), and hence, invariant under a change in the qubit's fiducial basis. Consequently, a 
judicious choice of the Schmidt basis gives D(tj(0)\ip),ij(l)\^})) = AjA 2 = (1 — trpj 2 )/2, where Aj and A 2 are the 
Schmidt coefficients in the decomposition between the j-th qubit and the remainder of the system. 

Brennen's simplification immediately allows the generalization of Q to multi-qudit states \ip) £ (C D )® n , and by 
considering all other possible bipartite divisions, we can now define a class of related multipartite entanglement 
measures in the obvious manner: 

^ /,n D m (, m\(n-m)\ ^ ,\ 

QmW = J^—J 1 , E tr ^ m = 1, . . ., \n/2\ (6) 

V "' \S\=m I 

where Sc {1, . . . , n} and ps = tr gi \ip}(tp\ is the density operator for the qudits S after tracing out the rest. Note that 
Q rn reduces to the original Meyer- Wallach measure when m = 1 and D = 2. The above "multipartite" entanglement 
measures are merely averages over the well-established bipartite measure. Consequently, Q m is invariant under 
local unitary transformations, nonincreasingon average under local quantum operations and classical communi- 
cation i.e. Q m is an entanglement monotone |43| . and < Q m < 1. The lower bound is only reached for product states. 

Proposition 1: Q m (V ; ) =0 iff \ip) = <8)j=i l^j) f° r some £ C D i.e. is a product state. 

When m = 1 the upper bound is reached by the generalized GHZ states 



M-^g lo- 



rn 
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In general 



>(7) = 1~ J)m _ 1 ■ (8) 



and hence, the entangled states I7) do not saturate the upper bound for m > 1. We have not, however, established 
whether or not there even exist states which saturate the upper bound. 

Define an m-uniform multi-qudit state to be a state with the property that after tracing out all but m qudits we 
are left with the maximally mixed state, for any m-tuple of qudits. Thus, all information about the system is lost 
upon the removal of n — m or more parties. 

Proposition 2: Qm(ip) = 1 iff ps = trg/ [V'KV'I = D~ rn l whenever |5| = m i.e. is m-uniform. 

Obviously, if \ip) is m-uniform then it is also (m — l)-uniform, and hence, Q m {ip) = 1 =>■ Q m -i(ip) = 1. However, 
note that the measures Q m do not obey any ordering. For example, in the case of qubits, consider the generalized 
W-states 



1 " 

^|q}«M-i ®|i)®|0)® n --*. (9) 



in 



One can calculate 



^ / x 2 m+1 (n-m)m , , 

and hence, for n = 6 say, Qi = 5/9 < Q3 = 4/7 < Q2 — 16/27. The measures Q m also do not preserve the 
partial ordering of entangled states i.e. Q m '(V0 < Q m '(4>) does not necessarily imply that Q m (i>) < Q m (4>) f° r 
other m. These facts might be considered as unlucky properties of Q m - However they do suggest that the extremal 
entanglement measure Q[ n /2] docs not necessarily tell the entire story; different Q m capture different aspects of 
multipartite entanglement. The original Meyer- Wallach measure Qi is the average entanglement between individual 
qudits and the rest, whereas, on increasing m, Q rn measures the average entanglement between blocks of qudits, of 
an increasing size, and the rest. Consequently as m increases, we expect that Q rn will be sensitive to correlations of 
an increasingly global nature. 

Proposition 2 implies that the task of finding states which saturate the the upper bound 1 of Q m is equivalent to 
the construction of m-uniform multi-qudit states. We now show in the next section how quantum-error-correcting 
codes (QECC's) produce m-uniform multi-qudit states. An example is the six-qubit hexacode state which arises 
as the code subspace of the self-dual qubit stabilizer code [[6,0,4]]. In this case Qi{H) = Q2{H) = Qz{H) = 1. 



III. MULTIPARTITE ENTANGLEMENT AND QECC'S 

The idea behind quantum error correction [22, 0, EH EE H3> EE E3 is to encode quantum states into qudits in 
such a way that a small number of errors affecting the individual qudits can be measured and corrected to perfectly 
restore the original encoded state. The encoding of a i^-dimcnsional quantum state into n qudits is simply a linear 
map from C K to a subspace Q of (C D )® n . The subspace itself is referred to as the code and is orientated in such a 
way that errors on the qudits move encoded states in a direction perpendicular to the code. 

A. General QECC's 

An error operator E is a linear operator acting on (C 15 )®". The error is said to be detectable by the quantum code 
Qif 

y>\E\il>) = (<fi\E\4>) (11) 

for all normalized l^), [</>) € Q. Equivalently, if Q is spanned by an orthonormal logical basis I 3 ' — 0, . . . , K — 1}, 
then an error E is detectable if and only if 

{j L \E\i L )=C{E)5ij (12) 
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for all < i, j < K — 1 where the constant C(E) depends only on E. It is a general theorem of QECC's that a set of 
errors £ can be corrected by a code Q, if and only if for each Ex,E% e £ , the error E\E\ is detectable by Q. 
A /oca? error operator has the form 

£ = Mi <g> • • • <g> M n (13) 

where each Mi acts on C 15 . The weight of a local error operator E, denoted by wt(E), is the number of elements Mi 
which are not scalar multiples of the identity. A quantum code Q has a minimum distance of at least d if and only if 
all local error operators of weight less than d are detectable by Q. A code with minimum distance d = 2t + 1 allows 
the correction of arbitrary errors affecting up to t qudits. In the case of qubits, such codes are denoted by the triple 
((n,K,d)). We will use the notation ((n 7 K,d))rj for the general case of qudits [5(|. An {{n,K,d))o code is called 
pure if (ip\E\ip) = D~ n trE for all S Q whenever wt(E) < d. When considering self-dual codes (K — 1), we adopt 
the convention that the notation ((n, l,d))u refers only to pure codes since the condition on the minimum distance 
is otherwise trivial. 

There is a continuum of possible errors in a single qudit; however, due to the phenomenon of measurement collapse, 
the correction of an arbitrary single-qudit error only requires an ability to correct D 2 different types, each correspond- 
ing to an orthonormal basis element for single-qudit operations. One choice for a nice error basis |5lL l58j | is the 
displacement operator basis 

D(/j,, v) = e "^/D x ^Z" < n, v < D - 1 (14) 
where the Weyl operators X and Z arc defined on a basis {\j) \ j = 0, . . . , D — 1} for C D through the equations 

X\j) = \j + lmodD), Z\ ] )=e 2 ^' D \ ] ). (15) 

The displacement operators reduce to the Pauli matrices for qubits, satisfy the relations 

D(n,v) = e l ™D{ii + D,v) = e^Dfav + D) (16) 

D{p,v)t = D(-n,-v) = ^+ V + D )D(D- h,D-v) (17) 

D(n,v)D(a,f3) = e 2mil/a -^ )/D D{a,P)D(n,is) = e wi ^ a -^ )/D D(fi + a, v + 0) (18) 

tr [D(ix,rfD(a,f3)] = DS^p, (19) 

and thus form an orthonormal basis for all single-qudit operators: 

D-l 

A=-Y,to[D(l*,v) i A]D( Ji ,v). (20) 

ii.v— 

Similarly, the operators 

T>(/j,, v) = V(m ...fin,ui... v n ) = D(ni, v{) (g • • • (g) D(fi n , v n ) < w,, Vk < D - 1 (21) 

form an orthonormal basis for the set of all n-qudit operators: A = D~ n v tr[2?(^t, vy A]T>(/j,, v). The weight of 
T>(fi, v) is simply the number of pairs v^) different from (0, 0). We are now in a position to make a more explicit 
definition of what we mean by an ((n, K, d))r> QECC. 

Definition: Let Q be a i^-dimensional subspace of (C D )®" spanned by the orthonormal logical basis {\]l) \ j = 
0, . . . , K — 1}. Then Q is called an ((n, K, d)) d quantum- error- correcting code if 

(j L \V( t i,u)\i L ) = C(n,u)6 ij (22) 

for all T)(n, u) with wt[I?(/x, i>)] < d and < i,j < K— 1. If C(/Lt, u) = S^oS^o the code is called pure. An ((n, 1, d))n 
code must be pure by convention. 

An ((n, K, d))n QECC can detect and recover all errors acting on < d/2 qudits. It is now evident how quantum 
codes produce maximally entangled states. 

Proposition 3: Q m (ip) = 1 iff \ip) is a (pure) ((n, 1, m + 1))b quantum-error-correcting code. 
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Proof: If Q m (ip) — 1 then is m- uniform, and consequently 

= tr[MMX>(/i,»')] (23) 
= D~ n tr [£>(/*,!>)] (whenever wt[T>(jj,,v)] < m) (24) 
= ^o^o (25) 

given that the displacement operators are traceless for all (/i, v) ^ (0, 0). Thus, \tp) is an ((n, 1, m + QECC. 
Conversely, if \ip) is an ((n, 1, m + 1))d QECC, then rewriting m the displacement operator basis 

D n \^\ = 1 + Y, ^V(n,v) + Y, c^Vfav) (26) 

l<wt[X!(/i,i/)]<m m+l<wt[X>(/i,i/)]<n 

we see that the coefficients — (ip\T>(fi,u)\ip) are nonzero only in the second sum, and hence, given the traceless 
property of the displacement operators, 

Ps = trg,\il>)M =D- m l (27) 
whenever \S\ = m. Thus \ip) is m-uniform and Q m {^) = 1. □ 

Note that any state G Q, where Q is a pure ((n, K,m + 1))d QECC, is itself an ((n, 1, m + 1))d QECC. 
Consequently, pure ((n,K,m + 1))d codes define entire subspaces of maximally entangled states. The connection 
between quantum codes and entanglement is noted in , and alluded to elsewhere 0>E(J; however, we cite the 
work of Rains |54| for a rigorous proof of the relationship even though no mention of entanglement can be found in 
the paper. Here quantum weight enumerators were studied extensively. It will later prove advantageous to now revisit 
Rains' work in the current article. 

Defining Pq as the projector onto the code subspace Q with dimension K, the Shor-Laflamme enumerators of a 
quantum code are (5^ 

MPq) = ^ E |tr[2?(/^)P Q ]| 2 (28) 

wt [D ([A 

= J? E tr[VOjL,u)PaD(p,u^P Q ] (29) 

wt \T> {^jl 

where i = 0, . . . , n. Rains [54[ denned two new enumerators 

4( p s) = ^2 E ^s[tr s/ [P Q } 2 ] (30) 

|S|=» 

B 'i( P Q) = Jl E tr *' [^s[Pq] 2 } (31) 

\S\=i 

related to the Shor-Laflamme enumerators via the equations 



i 'W 

A' m {P Q ) = D- m Yl w r, MPQ) (32) 

' (m — i)\(n — my. 

B' m (P Q ) = D- m Yl \u m V Bi ^- (33) 

■* — ' (m — i)\[n — m)\ 

i=0 v ' y ' 

This relationship was only given in the qubit case where the displacement operators reduce to Hcrmitian Pauli 
matrices. However the proof extends easily to qudits with the the help of Eq. I|17|) . It is easy to see that the weight 

enumerators satisfy the normalization condition A' q (Pq) = B' q (Pq) = A (Pq) — B (Pq) = 1, and for self-dual codes 
{K = 1) 

B[{Pq) = A'^Pq) [b^Pq) = MP Q ) ] (34) 
for all < i < n. In general, the weight enumerators satisfy |54| 

B[{Pq) > A'^Pq) > \Bi(P a )>MPQ)>0] ( 35 ) 



6 



for all < i < n. 



Theorem [54]: Let Q be a quantum code with associated projector Pq. Then Q has minimum distance of at least 
«/ iff 



B' d - X {PQ) = A' d _ x {P Q ) [Bi{P Q ) = MP Q ) for all < i < d 

and is pure iff 

B' d - X {PQ) = A'^Pq) 



D x - d n\ 



(d- l)\(n-d+ 1)! 



B t (P Q ) = A l (P Q ) = for all < i < d 



(36) 



(37) 



Proposition 3 is now immediately apparent since 
D" 



QmWO = 



D" 



1 



= 1 - 



D' 



m\(n — i)\ 
J n\(m — i)\ 



Ai{H){^\). 



(38) 



Noting that KA! i — B' n _ i , one can use the above theorem to derive bounds on the minimum distance for general 
quantum codes. In the case of ((n, 1, d))n QECC's the following conditions must hold: 



A'i — A' n 

A = 1 

Ai = 

Ai > 



< i < n 

< i < d 
d < i < n. 



(39) 
(40) 
(41) 
(42) 



When d — [n/2\ +1, these equations uniquely specify the weight distribution {A;}. Solving equations H39(l - (|41ll . we 
obtain 



Ai = 



E 



(-iy-i(D 2 i- n - 1) 



d < i < n 



and under the condition A^+i > 0, we find that we at least require 

n < 



2(D 2 — 1) if n is even 

2D(D + l)-l if n is odd 



(43) 



(44) 



for an ((n, 1, [n/2\ + 1))d QECC to exist. Consequently, Q|„/2J W < 1 f° r au IV') whenever Eq. il^ is not satisfied. 
For d < [n/2\ we must resort to linear programming techniques on equations (|39H - l|42f) to prove the nonexistence of 
((n, 1, d))n QECC's. However tighter bounds could be obtained by using the generalized quantum shadow enumerators 
|56ll57j |. We make no attempt at this task in the current article, but instead specialize to qubits where many examples 
of QECC's are already known. 



B. Stabilizer qubit QECC's 

An important class of quantum codes are the so-called additive or stabilizer codes [44ll45| . A stabilizer code is defined 
as a joint eigenspace of an Abelian subgroup S (called the stabilizer) of the error group 

/i/., Vk, A < D — 1}. When D is prime, these codes can be described by an (n — k) x n stabilizer matrix over GF(D 2 ) 
and are examples of {{n,D k ,d))n QECC's. The notation [[n, k, d]]n is then used, or simply [[n,k, d]] when D = 2. 

A classical additive code over GF(4) of length n is an additive subgroup C of GF(4) n . In the case of qubits, stabilizer 
codes correspond to classical additive codes over GF(4) 0]. This is shown as follows. Letting GF(4) = {0, l,u,cJ} 
where to — uj 2 = 1 + u>, we define the conjugate of a; € GF(4), denoted x, by the mapping = 0, 1 = 1, and ZJ = u>. 
Next define the trace map Tr : GF(A) ->■ GF(2) by Tr(af) = x + x 2 i.e. Tr(Q) = Tr(l) = and Tr(w) = Tr(ZU) = 1, 
and the trace inner product of two vectors x = x± . . . x n and y = yi . . . y n in GF(4) n as 



71 

x*y = ]TTr(^yl) € GF{2), 

i=i 



(45) 
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The weight wt(x) of x G GF(4) n is the number of nonzero components of x, and the minimum weight of a code C is the 
smallest weight of any nonzero codeword in C. Next, by defining the mapping $ : GF(A) n — > £ by ^(x) = 2?(</> _1 (x)) 
where </>(/i, v) = ujfi + Z3u, we can associate elements of GF(A) with Pauli matrices (uj — ► X, ZJ — * Z, 1 — > iXZ, 
— * J), addition of vectors over GF(4) n with multiplication of operators in £ (neglecting phases), and the trace inner 
product on GF(A) n with the commutator on £. 

If C is an additive code, its dual is the additive code C = {x G GF(4) n |x*c = OVc G C|. The code C is 
called self- orthogonal if C C C and self-dual if C = C . The following theorem now applies |44| |: Suppose C is a 
self-orthogonal additive subgroup of GF(4) n , containing 2 n ~ fc vectors, such that there are no vectors of weight < d 
in C ± \C. Then any joint eigenspace of <&(C) is an [[n, k, d]] QECC. 

We say that C is pure if there are no nonzero vectors of weight < d in C x . The associated QECC is then pure if 
and only if C is pure. By convention, an [[n, 0, d]] QECC corresponds to a self-dual additive code C with minimum 
weight d. Consequently, [[n, 0, d]] QECC's are always pure, and are examples of ((n, 1, d)) QECC's which saturate the 
entanglement measures Q m . 

The advantage of making the above correspondence is that a wealth of classical coding theory immediately becomes 
available. Indeed the classical self-dual additive hexacode with generator matrix 



1 1 uj uj 

1 uj 1 uj 

1 uj uj 1 

u> uj ZJ ZJ 

lo ZU uj uj 

uj ZJ ZJ uj 



(46) 



gives the quantum hexacode [[6,0,4]] mentioned previously. The rows of the generator matrix define a basis (under 
addition) for the classical code C, and, with the above correspondence, define generators (up to a phase) for the 
stabilizer S in the quantum version. Another example is the [[2, 0, 2]] qubit code generated by 



1 1 

UJ UJ 



In this case the quantum code is an EPR state e.g. (|00) + |11))/V2. 
deleting the first row and column of the hexacode generator matrix [Eq. 
A [[3, 0, 2]] code 



(47) 



We can obtain an [[5,0,3]] qubit code by 
PBjl]- This process is called shortening |60| . 



1 1 

UJ UJ UJ 

1 1 



(48) 



is obtained by lengthening the [[2, 0, 2]] code. 

The four codes of lengths n = 2, 3, 5 and 6 mentioned thus far all produce quantum stabilizer codes with the 
property Q\ n /2\{'4>) — E Unfortunately, known bounds on such codes prevent this from being the case for other 
lengths. An additive self-dual code is called Type II if all codewords have even weight, and Type / otherwise. It can 
be shown that all Type II codes have even length. If dj, du is the minimum weight of an additive self-dual Type I, 
Type II code, respectively, of length n > 1, then |58ll59ll60| 



d 



ii 



< 



< 



2[n/Q\ + 1 
2[n/6j +3 
2[n/6j +2 

2[n/6j +2. 



if n = mod 6 
if n = 5 mod 6 
otherwise 



(49) 
(50) 



If a code meets the appropriate bound it is called extremal. A code is called optimal when it is not extremal and no 
code can exist with a larger minimum weight. The above bounds imply that [[n, 0, / 2j + 1]] stabilizer codes may 
exist only when n — 2, 3, 5, 6 and 7. However [[7, 0, 3]] codes are known to be optimal and we a left with the remaining 
four cases. 

The weight distribution of an additive code C 



A, 



|{x G C | wt(x) = i}\ 



(51) 



is also the weight distribution for the corresponding quantum stabilizer code, and thus, the entanglement of the 
stabilized state is easily calculated through formula (|38|l . We can see this by noting that for stabilizer codes, the 



8 



n 


(/ 
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Ai 


A 2 


A 3 


A 4 


^5 


A 6 
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Qz 
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2 


8 
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1 
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5 


3 


1 








10 


15 
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1 
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6 


4 
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1 


1 






7 


3 


1 








7 


21 


42 


42 


15 














1 


1 


34/35 






7 


3 


1 








3 


29 


42 


34 


19 














1 


1 


242/245 






8 


4 


1 











42 





168 





45 












1 


1 


1 


24/25 




8 


4 


1 











26 


64 


72 


64 


29 












1 


1 


1 


512/525 




9 


4 


1 











26 


48 


136 


160 


93 


48 










1 


1 


1 


932/945 




9 


4 


1 











18 


72 


120 


144 


117 


40 










1 


1 


1 


104/105 




10 


4 


1 











30 





300 





585 





108 








1 


1 


1 


104/105 


212/217 


11 


5 


1 














66 


198 


330 


495 


550 


330 


78 






1 


1 


1 


1 


216/217 


12 


6 


1 

















396 





1485 





1980 





234 




1 


1 


1 


1 


1 146/147 


13 


5 


1 














15 


236 


356 


1197 


1530 


2012 


1956 


650 


239 


1 


1 


1 


1 


13294/13299 26938/27027 



TABLE I: The weight distributions Ai and corresponding entanglement Q m for extremal (or optimal for n = 7 and 13) additive 
self-dual codes. In all but the cases n — 10 and 13 these are the only possible weight distributions. 



projection onto Q is given by 

P 2 = ||| E H e ) 1e (52) 

where S is the stabilizer and X(E) is the eigenvalue associated with E i.e. E\ip) = X(E)\ip) for all S Q. We remark 
that the quantum weight distribution Bi corresponds to the classical weight distribution of the dual code C^~. 

In Tabled the weight distributions for extremal (or optimal for n = 7 and 13) additive self-dual codes is collected 
I n a H but the cases n = 10 and 13 these are the only possible weight distributions. The weight distribution 
is unique for extremal Type II codes. Also tabulated is the corresponding entanglement Q m for the quantum code. 

Although the bounds mentioned above [Eq.'s Q49l50[l ] were given in the context of stabilizer codes, they also apply to 
general QECC's |^. Consequently for qubits, there exist states with <5[n/2j (VO = 1 on ly m the cases n — 2,3, 5, 6, 
and possibly when n = 7 where a nonadditive ((7, 1,4)) code might still exist. It is not known what the suprcmum 
of Q m {tp) is in general; however, given the examples in Tabled we expect it to be very close to 1 when n is large. 
It is interesting that the mean of Q m (i(>) (given in the next section) seems unaffected by the erratic behavior in the 
supremum. 

For the most part, quantum coding theorists have primarily studied qubit codes. Some work on qudit codes exists 
[50L Ir33l IrJil I60I Irjrj ]. but there are very few known examples. One exception is the generalization of the hexacode. 
A ((6, 1,4)).d code is known to exist for all D [Hfl. The ((6, 1,4)) d code belongs to a class of optimal codes called 
maximum distance separable (MDS) codes which saturate the quantum Singleton bound [5(j i-e. K = ]j n - 2d + 2 . 
Quantum MDS codes must be pure, and thus, define subspaces of maximally entangled states. Self-dual quantum 
MDS codes have weight distributions specified by Eq. J33J). Other examples of MDS codes include the [[6, 2, 2]]n and 
[[7, 3, 3)]d stabilizer codes which exist for all prime D 65]. More recently, the existence of some families of quantum 
MDS codes was proven [6rl |. For example, when D is a prime power and n is even, a self-dual ((n, 1, n/2 + 1))d code 
exists for all 3 < n < D. 



IV. MULTIPARTITE ENTANGLEMENT OVER SUBSPACES 

Using Lubkin's formula j63| for the average subsystem purity, one can easily calculate the mean entanglement for 
random pure states sampled according to the unitarily invariant Haar measure d/j, (J d/i(ip) = 1): 

f D m + 1 

(QmW)^ = / dM)Q m ty) = 1 - -^j- (53) 
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This shows that when the overall dimension D n is large, a typical state has nearly maximal entanglement. One could 
also consider the average entanglement over a subspace V determined by the projector Py: 



(Qm(i>)) i}ev = I dUvMQmM- ( 54 ) 



Proposition 4: 



where if = trPy = dimV. 



Proof: Consider an arbitrary bipartite system Tt = C Da ® C^ 3 and define the swap operators (1 < i < j < 4) 
which transpose the i-th and j-th factors of H® 2 . Using the identity tr[(A ® P)T] = tr[AB], where T is the swap, we 
first rewrite the subsystem purity of a state \tp) € W as 



trp A 2 = tr[|^)(V>|® 2 r 13 ] (56) 
where = tr^ l^*)^. Now consider the operator 

w = / ^vWI^^I 02 (57) 



supported on the totally symmetric subspace PV® 2 , where the projector P = (1 + Ti 3 T24)/2. If we choose ci/xy to be 
the unitarily invariant Haar measure on V, then [£/® 2 ,u] = for all unitary operators U E V(DaDb)- And since the 
group elements U® 2 act irreducibly on PV® 2 , by Schur's lemma u> is simply a scalar multiple of the identity (on 
PV m ). Hence, 



P v m P (58) 



+ 1) 

on H® 2 , where the constant factor is found through the normalization condition trw = 1. Thus 

J^dM4>)^PA 2 = R J + ^ tr [p v 02 (i + r 13 T 24 )r 13 ] (59) 

= [fl^is] +tr [Py« 2 T 24 ]) (60) 

1 



+ 1) 



trpV + trpV) (61) 



where /5 a = tr# Py and ps = tr^Py. We have derived the average purity over a subspace for an arbitrary bipartite 
system. Given that the measures Q m are simply averages over bipartite purities, one can now deduce the final result. □ 

In particular, for a QECC Q, one can now explicitly determine the average entanglement of encoded states in terms 
of the weight distribution of the code. For example 

= j^—^l - n ^ K ^i KA U P Q )+B: n (P Q )}) (62) 



= 1 - 



1 m\(n — i)\ 



(D m - 1)(K + 1) n\(m - i)\ 



^ it KAi (P Q ) + B t (P Q ) . (63) 



Such formulae make explicit the importance of entanglement as a resource for quantum error correction. Pure codes 
(Bi — Ai — 0, < i < d) necessarily have high levels of entanglement, but impure codes (Bi — Ai > 0, < i < d) 
need not. However the most compact codes (least n) all seem to be pure |55|. The relationship between entanglement 
and QECC's remains relatively unexplored in the literature. We do not, however, pursue this line of research any 
further in the current article. Elements of the proof of Proposition 4 were borrowed from the work of Zanardi et al. 
|23j where the concept of entangling power was defined. In the next section we investigate multipartite entangling 
power with respect to the measures Q m - 
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MULTIPARTITE ENTANGLING POWER 



Following the work of Zanardi et al. '23 |. we define the multipartite entangling power of the unitary operator 
U 6 U(D n ) acting on (C D )® n as simply the average entanglement generated over all product states: 

e p (U) = J dpntyu ■ ■ ■ , i> n )E(U\ih) ® • • • ® |Vn)) (64) 

where |^) € C D . The measure is chosen to be the product of n independent Haar measures over the 
constituent subsystems C D . Consequently, the entangling power is invariant under the action of local unitaries: 
e p (Ui (8 • • • (g) U n UV\ ® • • • ® 14) = e p (/7) for all [/,, V$ € U(-D). If we now restrict our attention to the entanglement 
measures E = Q m , the calculation of e p is facilitated by a simple formula. 



Proposition 5: 



e^(U) = - 



D' 



1 - 



m\(n — to)! 



(65) 



where the average subsystem purities 



Rs(U) = 



D(D + 1) 



tr 



|S|=m 



(66) 



\i=l / Vies 

the swap operators Ty (1 < i < j < 2n) transpose the i-th and j-th factors of (C D )® 2n , and Pjj = (1 + Py)/2. 



Proof: The derivation of this formula is similar to that for Proposition 4 and follows Zanardi et al. 23]. We first 
rewrite the subsystem purity of a state |\&) £ (C D )® n as 



By choosing |<f) = U\ipi) 



tr ps 2 = tr (|¥><* | ® 2 l[T ul+ \ 

\ip n ) we have 

R S {U) ee dn n (ipi,...,ip n )trps 2 



trf^ 2 fiC7t®2-Q T .. + J 



where 



O = / d(in(l/>i, ■ ■ ■,1pn)(\ll>l){'>l>l \ ® • • • ® |V>n)(V>n 



Z'2 



(67) 



(68) 



(69) 



(70) 



Now, considering the operator w = j d^ii^i)]^) (ip\® 2 supported on the totally symmetric subspace Pi2(C D )® 2 , where 
Pi 2 = (l + T 12 )/2, we know from previous results (see proof of Prop. 4) that u = 2/D(D + l)P 12 on (C D )® 2 . Finally, 
given that Q factorizes into the product of n independent averages of the form ui, we have 



= 



D(D 



it) n* 



-i+n 



(71) 



and our final result. □ 

Given our definition of the entangling power [Eq. I|64(l ]. the (Haar measure) average of e p (U) over U(D n ) is 
equivalent to the average entanglement found in random states: 



1 



1 



Thus, typical unitaries generate nearly maximal entanglement when the overall dimension D n is large. 



(72) 
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FIG. 1: (a) Phase-space portraits of the classical kicked rotor for k — 0,0.2,1 and 6 (from top to bottom), and (b) the 
corresponding entangling power in the quantum case. A Hilbert space of 6 qubits was chosen allowing the investigation of the 
multipartite entanglement measures Q m for m = 1 (dotted), 2 (dashed) and 3 (full). The average entanglement for random 
states is shown in the lighter tones. 



VI. AN APPLICATION AND CONCLUSION 



imediate ap plication ot Proposition 5 (and 1231) occurs in the study of the entangling capabilities of chaotic 
systems H3, HE Ea 123, [H, [2^ E3, EH H3, HE [33, Hall HHH HI - Consider a classical map of the toroidal phase 
space [0, l) 2 . A quantized version may be constructed in a Hilbert space of dimension N spanned by the position 
states \qj), where qj = (j + 1/2) /N and j = 0, . . . , N — 1. By choosing N = D n , we can map our Hilbert space onto 
the tensor-product space (C D )® n through the correspondence 



Fra) ' 



i=l 



€{0,. ..,£>-!}, 



(73) 



and hence, use the measures Q m to investigate the quantum map's multipartite entangling power. The different 
constituent qudits Xi address the coarse (small i) and fine (large i) scales of position. Consequently, for chaotic 
maps where phenomena such as mixing and exponential sensitivity are generic, we expect high levels of entanglement 
generation. This was noted in [36] where the entangling power of the quantum baker's map |69l ItoL l7ll l72l l73l | was 
investigated. 
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For example, consider the kicked rotor (or standard map) 74] : 

q n +i = Qn+Pn+i modi (74) 
k 

Pn+i = Pn + — sin 2nq n modi. (75) 
The entangling power of the quantum version 

N-l 

U\q 3 ) = e _i ^ cos27r9j ^ e^ ( '- j)2 |g;) (iV even) (76) 
;=o 

constructed in a Hilbert space of 6 qubits (N — 2 6 ) is plotted in Fig. ^b). Here we choose the parameter values 
k = 0, 0.2, 1 and 6 (from top to bottom), corresponding to the classical phase spaces drawn in Fig. QJa). As expected, 
the entanglement saturates at a value predicted for random states [Eq. (|53[) ] upon the appearance of chaos in the 
classical map. An alternative interpretation may be that quantized chaotic maps produce unitaries whose powers are 
typical in the space of all unitaries. This follows from Eq. I|72l) . 

In conclusion, we have shown that the average bipartite entanglement, Q m , is a useful measure of multipartite 
entanglement, presenting a relationship between these measures and quantum-error-correcting codes. This was done 
by deriving an explicit formula relating the weight distribution of the code to the average entanglement of encoded 
states. We have also extended the work of Zanardi et al. |2^| on entangling power to the multipartite case. Although 
the entanglement measures considered in this paper provide little intellectual gratification, their simplicity allows 
perhaps more important attributes such as computability and applicability. We must stress, however, that in defining 
such simple measures we offer no progress towards a deeper understanding of the nature of entanglement in the 
multipartite case. 
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